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ABSTRACT 

We present the luminosity function of [OII]-emitting galaxies at a median 
redshift of z = 0.9, as measured in the deep spectroscopic data in the STIS 
Parallel Survey (SPS). The luminosity function shows strong evolution from the 
local value, as expected. By using random lines of sight, the SPS measurement 
complements previous deep single field studies. We calculate the density of in- 
ferred star formation at this redshift by converting from [Oil] to Ha line flux as 
a function of absolute magnitude and find p = 0.043 ± 0.014 M Q yr _1 Mpc~ 3 at 
a median redshift z ~ 0.9 within the range 0.46 < z < 1.415 (H = 70 km s _1 
Mpc" 1 , Qm = 0.3, = 0.7). This density is consistent with a (1 + z) 4 evolu- 
tion in global star formation since z ~ 1. To reconcile the density with similar 
measurements made by surveys targeting Ha may require substantial extinction 
correction. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
fundamental parameters — 
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1. Introduction 

The comoving density of global star formation in the Universe has decreased significantly 
since a redshift of z ~ 1 — 2 (Hogg 2002; Madau, Pozzetti, & Dickinson 1998; Madau et 
al. 1996), but the details of the evolution are still uncertain. Extensive measurements of 
the star-formation history have concentrated on surveys for Ha emission (e.g. Gallego et 
al. 1995, Yan et al. 1999, Tresse et al. 2002) and the UV continuum (e.g. Lilly et al. 1996, 
Connolly et al. 1997). A correction factor is needed reconcile the two methods (Calzetti et al. 
1997; Tresse & Maddox 1998; Glazebrook et al. 1999) as dust obscuration is expected to hide 
some significant fraction of star formation in the distant Universe. Complementary optical 
data is provided by targeting the [Oil] A3727 emission line, which can also provide a direct 
measurement of the rate of star-formation, though its calibration depends on the metallicity 
and ionization state of the gas, as well as the reddening (Kennicutt 1998). Indeed, surveys 
at higher redshift will need to depend on the oxygen line. At z > 4, the H/3 line redshifts 
out of the K-band, making it inaccessible to ground-based NIR spectrometers. As a result, 
[Oil] emission is the best rest-frame optical tracer of star-formation available at high redshift 
until NGST. 

In this paper, we present analysis of the deepest spectroscopic fields in the STIS parallel 
survey (SPS; Gardner et al. 1998). Our motivation is to determine the comoving density of 
star- format ion at z ~ 0.9 as inferred from the [Oil] emission line. The deep SPS data cover 
141 square arcminutes, and detect objects down to continuum magnitudes of m^occD ~ 26 
in filterless imaging, with emission line strengths as faint as 7.5 x 10~ 17 ergs cm' 2 s _1 
(Teplitz et al. 2002, hereafter Paper I). These data are complementary to other surveys of 
the [Oil] luminosity density. Hogg et al. (1999a; hereafter H99) measured [Oil] emission in 
375 galaxies with R < 23.3 from the redshift survey data in the 50 square arcminute "Caltech 
area" (Cohen et al. 2000) around the Hubble Deep Field North (HDF; Williams et al. 1996). 
The HDF spectra obtained with LRIS (Oke et al. 1995) on the Keck telescopes measure 
emission lines to a limit 15 times fainter than the SPS. Prior to H99, the best measurement 
of the [Oil] luminosity density was obtained the Canada France Redshift Survey (CFRS; 
Lilly et al. 1995) by Hammer et al. (1997, hereafter CFRS9). A significant increase in [Oil] 
luminosity is seen with lookback time in the range < z < 1 (Cowie et al. 1996; Gallego et 
al. 2002). 

Throughout the paper, we will adopt a flat, A-dominated universe (H = 70 km s -1 
Mpc" 1 , VL M = 0.3, = 0.7), except where otherwise noted. 
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2. The Data 

The slitless spectroscopic capability of STIS, using the G750L grating with spectral 
resolution of ~ 750 at 7500A, allowed a random-field redshift survey for star-forming galaxies 
at 0.46 < z < 1.7. The deepest 219 of the ~ 2500 fields, which were observed by STIS in 
parallel with targeted HST observations, contain 78 objects identified with [Oil] emission 
(Paper I). Of these, 14 have multiple emission lines, providing secure redshifts. For the 
rest, the identification of a single line as [Oil] is based on the non-detection of other lines, 
together with photometric information from the paired direct filterless imaging. A signal to 
noise limit of a — 2.5 was imposed for the strongest line in an object in order for it to be 
included in the catalog. For the present analysis, we impose a redshift cutoff of z < 1.415 to 
avoid the objects in the poor sensitivity region at X b s > 9000 A. In total, we include 71 of 
the 78 [OII]-emitters from Paper I. 

The SPS detections are biased towards objects with a strong contrast between line and 
continuum. This bias means most SPS line-emitters are morphologically compact and that 
most of the emission lines have high equivalent width (EW). The SPS detection limit is 
~ 35 A for observed EW, consistent with the values for most galaxies undergoing rapid 
star-formation (Kennicutt 1992). The measured EW value is uncertain for most SPS objects 
due to low significance of continuum measurement in the spectra. Typically the continua are 
measured with a signal-to-noise ratio (SNR) of a c ~ 0.5 per two-pixel resolution element. 
The continuum is measured over many resolution elements (typically 20 — 30), but is still a 
major source of error in the EW. Figure 1 compares the EW distribution for the SPS with 
that of other [Oil] surveys at similar redshifts. The SPS objects have a median rest-frame 
EW for [Oil] of W (OII) ~ 70 A, but many of the W (OII) > 100 A objects have poor 
continuum detections, so that the EW may be overestimated (see H99). It appears from 
the figure that the SPS is missing a substantial portion of low EW lines, particularly at low 
redshift. At high redshift, low EW lines may simply have been missed in all of the surveys. 
Nonetheless, CFRS9 and Cowie et al. (1996) both conclude that there is a strong increase 
in the fraction of galaxies with strong [Oil] emission (EW([OII}) > 25 A) at z > 0.5. 



2.1. Incompleteness 

The SPS survey detects W(OII) > 35 A in the observational frame. Figure 2 shows 
the distribution of missing low-EW objects as a function of [Oil] luminosity, based on H99. 
In the analysis of the SPS sample, we correct for this incompleteness. Because the H99 
published data only present the equivalent width distribution for objects with continuum 
detection (2a), this incompleteness correction is strictly an upper limit. 
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The primary source of incompleteness in the SPS sample is missing objects that would 
have met the selection criteria. Visual detection by multiple observers was more accurate 
than any available automated detection scheme, but human error remains a factor. Simulated 
spectroscopic images were used to determine the completeness level of our detection scheme. 
Test images were created by inserting known spectra into SPS frames and trying to recover 
them. Two spectra were used as the basis of this experiment, one with W(OII) ~ 100 and 
the other with W(OII) ~ 35. Both spectra were high signal to noise (a n > 8), but were 
scaled to lower significance in the test frames. The emission line itself was moved within 
the spectrum to prevent identification by pattern matching. The spectra were inserted at 
random spots in randomly chosen images and scaled to an emission line signal-to-noise ratio 
of 3-7. For each test case, 300 simulated frames were created. In some frames, continuum 
spectra were also added without an emission line to further disguise the simulated data. 
Table 1 lists the resulting incompleteness as a function of EW and SNR. The completeness 
of the survey at any given SNR and EW value is obtained by interpolating between the data 
points in the table. 



3. Calculating the Luminosity Function 

We use the 1/Vmax method to calculate the luminosity function (LF) of [Oil] emission 
lines. For each galaxy we determine the comoving volume over which it could be detected 
in our survey. The volume is different for each of the parallel fields, based on the limiting 
flux. Assuming we have j detections in k fields, then for a single galaxy in a single field we 
have, following Hogg (1999b) : 

W,,* =(^r) f r X (1 t? } ° D W (i) 



.H J J J z E(z) 

where, 



e{z) = y/n M (i + z)3 + n k (i + z)* + n A (2) 

and, Da is the angular diameter distance, which in the case of a flat, matter dominated 
universe is: 

/ c \ 1 f z dz' 

DA -ywJoTz)J EW) (3) 

Further, dQ is a function of z, because the available area depends on the placement of 
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the observed wavelength of the line and the position of the object within the field of view. 
That is, if the object is too far to one side, the line will fall off the edge of the detector. So 
for an npix^, x npix y detector, a central wavelength A c , a line at rest wavelength A and a 
dispersion d, 



f > n / • \X (1 + z) - \ c \\ . 

/ dfl = [npix x Jnpix^ (4) 

The effective area of the detector is slightly smaller than the default field of view 
(51.2"x52.2"), since we only consider the region of the registered frame with the full ex- 
posure time. Similarly, we do not consider area blocked by very bright stars or local galaxies 
(see the first table in Paper I). 

The volume integral is over z m i n to z max , as defined by the signal-to-noise in the line. 
In an exposure with limiting flux fu m , a line that is detected at X b. s with flux f b s and 
originating at z obs would be have been observable at any z such that 

where S(z) is the sensitivity function of the spectrum, normalized to unity at \ Q b s , and D L 
is the luminosity distance 



The individual Vmaxj,fe for the each of the survey fields are summed for each galaxy 
to obtain its Vmaxj- Finally, the number density of galaxies in a luminosity bin of width 
A(logL) is obtained by: 



where 77 is the inverse of the completeness function at the SNR and EW of the detected 
line. Figure 3 shows the SPS [Oil] luminosity function compared to that of H99 and the 
local measurement (Gallego et al. 2002). The SPS measurement appears systematically 
underdense at high luminosities. Only two objects contribute to the last log L Q n bin, so 
missing or misidentified objects may be the cause of that bin being so extreme an outlier. 
The flattening of the H99 LF at low luminosities compared to the SPS measurement might 
be expected given the H99 magnitude limit, since H99 do not detect the high-EW R > 23 
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galaxies that make up much of the SPS sample and may contribute to the faint end of the LF. 
The LF inferred from SPS shows strong evolution with respect to the local LF, as expected. 

4. Comoving Density of Star-Formation 

Ha emission is a good indicator of star formation because it traces the ionizing flux 
from hot stars. To infer the star formation rate (SFR) in a galaxy from [Oil] emission, it is 
necessary to assume an [Oil]: Ha ratio. This ratio has an average of 0.45 for local galaxies 
(Kennicutt 1998), but is highly dependent on the metallicity and reddening of the individual 
galaxy. Jansen, Franx, & Fabricant (2001) show that the ratio can vary by up to a factor 
of 7, and that it has a strong inverse correlation with continuum luminosity. We adopt the 
relation they find for local galaxies with strong Ha emission (EW(Ha)> 10 A): 

log(I^l) = 0.09M B + 1.42 (8) 
Ha 

Jansen et al. find considerable scatter about this relation. Apparently, the relation holds at 
higher redshift (z ~ 1; Tresse et al. 2002), but it has not been tested on many objects, and it 
may not hold in all cases (Hicks et al. 2002). If we instead use the average value for the ratio 
of [Oil] to Ha, then we obtain similar results though the distribution of Ha luminosities has 
less outliers; that is, we find fewer galaxies with L^ a < 10 42 or > 10 43 . 

In order to apply the Jansen et al. relation, it is necessary to know the rest-frame Mb 
of each galaxy. The SPS images are taken with the STIS filterless CCD, which admits a 
much wider wavelength range than the B filter. We adopt a continuum slope proportional 
to wavelength, except in cases where a good detection of the continuum in the spectrum 
allows direct measurement. This conversion is an important systematic error in the SPS 
data, and in the future could be improved with photometry of the emission line objects in 
several filters. The average value of [OII]:Ha remains close to 0.45 for the SPS sample using 
the Jansen et al. relation, but some individual galaxies change by a factor of more than two. 

Figure 4 shows the luminosity function of SPS emission lines converted to Ha, compared 
to other measurements of LF(Ha) at similar redshifts, from the literature. Of particular 
interest are the comparisons to measurements made with NICMOS slitless spectroscopy on 
HST (Thompson et al. 1998). The NICMOS parallels (McCarthy et al. 1999; Yan et al. 
1999) are a survey comparable in many respects to the SPS, since they are biased towards 
compact high-EW objects and sensitive to a similar range of luminosities. The NICMOS 
parallels do not detect lines with EW Tes ^(Ra) < 50A, but Ha EW is typically 2.5 times that 
of [Oil] (Kennicutt 1998). 
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The SPS measurement of the Ha LF appears systematically somewhat lower than that 
of NICMOS. The SPS objects are at a slightly lower average redshift, and some evolution 
in the LF is expected. In addition, [Oil] is subject to more extinction than Ha, though 
this effect is mitigated by reddening in the galaxies used to fit the Jansen et al. relation for 
[OII]:Ha. 

We fit a Schechter function to our LF(Ha), even though the SPS fields do not contain 
enough detections to constrain all three parameters (<f)*,L*, and a) accurately. We adopt the 
faint-end slope measured for the local universe (a ~ —1.35; Gallego et al. 1995), as do most 
of the Ha surveys. Hopkins, Connolly, & Szalay (2000) use a steeper slope, a = —1.6, in 
good agreement with that of Sullivan et al. (2000) for [Oil]. If we had used this value, our 
final SFR densities would have been ~ 25% higher. The integral of the Schechter function 
is given analytically by 

0*LT(2 + a). (9) 

Thus the integrated Ha luminosity density inferred from the SPS data is 5.9 x 10 39 /i7 ergs 
s _1 at a median redshift of z ~ 0.9. This value is about 4.1 times the local Ha luminosity 
density (Gallego et al. 1995), but only half that measured in the NICMOS parallels at z ~ 1.3 
(Yan et al. 1999). 

Assuming case B recombination, and a Salpeter initial mass function truncated at 0.1 
and 100 M , the star formation rate in a galaxy can be inferred from the Ha luminosity 
using the Kennicutt (1998) relation: 

SFR(M G yr- 1 ) = 7.9 x 10~ 42 L(Ha) (ergs s" 1 ). (10) 

Using this conversion, we obtain a comoving density of star formation of 0.043 ± 0.014 
M y^ 1 Mpc~ 3 at z ~ 0.9. Figure 5 compares the SFR density to values measured from 
other rest-frame UV-optical surveys. The SPS value is calculated as an integral over a large 
redshift range (0.46 < z < 1.415), but it can also be measured for a subset of the sample. In 
Figure 5, we show the SPS density of star formation for the subsets of the data in the ranges 
0.46 < z < 1 and 1 < z < 1.41, with median redshifts of 0.72 and 1.14, respectively. We 
do not plot the SFR density inferred by H99, but these authors note that their SFR density 
agrees with CFRS9. 



5. Discussion 

As expected, the SPS data show strong evolution in the density of star-formation over 
the range 1.14 > z > 0.72. Although the uncertainties are involved are large, the evolution 
observed is consistent with (1 + z) 4 , the same as that computed by Tresse et al. (2002) for 
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evolution in Ha. The same exponent is derived when the other [Oil] measurements in Figure 
5 are included. The value agrees with the wavelength-independent exponent of 3.8 ± 0.8 for 
a matter dominated flat Universe (Hogg 2002). However, Baldrey et al. (2002) set a strong 
upper limit on the exponent of < 2 at z > 1. 

It is apparent from Figure 5 that the [Oil] measurements of the SFR density do not 
agree with the Ha measurements any more than with those inferred from the UV contin- 
uum. The clear difference between the Ha and UV measurements is usually attributed to 
dust extinction (e.g. Yan et al. 1999), but at 3727 A, the [Oil] emission should be subject to 
less extinction from dust than the UV continuum. Furthermore, any difference in extinction 
between Ha and [Oil] should have been included, to first order, in the calibration of the 
ratio of the lines by Jansen et al. (2001), as that ratio was derived empirically from typically 
metal rich galaxies. If more sensitive spectro-photometry were obtained for the SPS sample, 
it would be possible to recalculate the inferred Ha luminosities using reddening information 
in the Jansen et al. fit. It might be argued that the SPS selection effects favor more red- 
dened objects, if reddening is greater in more compact starbursts, perhaps as a result of the 
greater concentration of gas and dust. The CFRS9 measurement, however, shows the same 
distinction, and that survey had very different selection criteria. 

The disagreement between the [Oil] and Ha values might be resolved if both are cor- 
rected for extinction. Gallego et al. (2002) find the extinction correction for [Oil] to be 
a factor of ten on average, which brings their measurement of the local SFR density into 
agreement with the density inferred from Ha (Gallego et al. 1995). This value seems quite 
large, however, and does not agree with that used in other studies. Tresse et al. (2002) apply 
a factor of 2 extinction correction to Ha and only a factor of 3.6 to [Oil], similar to Sullivan 
et al. (2000). 

The uncertainties in the measurement of the SFR density are still large, and the number 
of objects with both [Oil] and Ha line fluxes is still small. A direct approach to resolving the 
differences suggested by Figure 5 is additional multiwavelength study of consistent samples. 
At z ~ 1, these measurements require sensitive near-IR spectroscopy. Observations of that 
kind are possible and becoming more common (e.g. Pettini et al. 2001, Teplitz et al. 2000) 
from the ground, or soon with the refurbished NICMOS. 
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Table 1. Incompleteness in the SPS 



SNR (W obs ~ 100) (W obs ~ 35) 

Completeness (%) Completeness (%) 

3 39 20 

4 56 28 

5 72 44 

6 86 58 
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Fig. 1. — The redshift distribution of rest-frame [Oil] equivalent widths in the SPS. The 
small filled circles are SPS data with relatively poor continuum spectra (SNR < 1 in two 
pixels), and the large filled circles have good continuum detections (SNR > 1 in two pixels). 
The measurements in the HDF (open squares; H99) and the CFRS (open triangles; CFRS9) 
are shown for comparison. The dashed line indicates the W Q b s > 35 A cutoff of the SPS 
data. 
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Fig. 2.— The histogram of [Oil] luminosities in the HDF from H99 (H = 100 km s" 1 
Mpc -1 , Vt M = 0.3, JIa = 0). The solid line is the entire sample, and the dashed line is 
the fraction of the sample with observed EW> 35A; i.e., those EW greater than the SPS 
equivalent width cutoff. 
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Fig. 3. — The luminosity function of [Oil] emission in the SPS compared to the 0.35 < z < 
1.5 sample in the HDF (H99). The open circles are the raw SPS detections. Small filled 
circles have the incompleteness correction applied down to the EW cutoff and the large filled 
circles have the additional correction for missing objects below the EW limit (see text). The 
open triangles are the [Oil] LF for the local Universe (Gallego et al. 2002) and the open 
squares are the HDF (from H99). 
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Fig. 4. — The luminosity function for [Oil] line fluxes converted to Ha (see text) compared to 
other measurements at similar redshifts. The large filled circles are the SPS values corrected 
for incompleteness including the missing low EW objects. The other data points are as 
follows - open diamonds: z ~ 1.3 from the NICMOS parallels(Yan et al. 1999); open squares: 
z ~ 1.25 from NICMOS observations of the Groth strip (Hopkins et al. 2000), each bin 
plotted for both spectroscopically confirmed (lower) and all candidate (upper) objects; open 
upward triangles: z ~ 0.73 from the CFRS (Tresse et al. 2002). The solid line shows the 
Schechter function fit with a = —1.35. 
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Fig. 5. — Comoving density of star formation as a function of redshift. Gray filled squares 
show the values derived from the Ha measurements (in order of increasing redshift) of: 
Gallego et al. (1995); Tresse & Maddox (1998); Pascual et al. (2001); Tresse et al. (2002); 
Glazebrook et al. (1999); Hopkins et al. (2000); Yan et al. (1999). The open squares are 
SFR densities inferred from UV continuum measurements (without extinction correction) 
by: Treyer et al. (1998); Lilly et al. (1996); Connolly et al. (1997). The small filled circles 
are the [Oil] luminosity density of Hammer et al. (1997; converted to Ha by [OII]/Ho;=0.45, 
see Kennicutt 1998) and Sullivan et al. (2000). The large filled circles are the integrated 
SPS luminosity function. The three points are not independent as indicated by the redshift 
ranges (horizontal bars) - the middle point includes all objects in the survey, while the other 
points are for the z < 1 and z > 1 subsamples. 



